Beating of Friedel oscillations induced by spin-orbit interaction 



S. M. Badalyan, 1 - 2 ^ A. Matos-Abiague, 1 G. Vignale, 3 and J. Fabian 1 

department of Physics, University of Regensburg, 93040 Regensburg, Germany 
J Department of Radiophysics, Yerevan State University, 1 A. Manoukian St., Yerevan, 375025 Armenia 
'Department of Physics and Astronomy, University of Missouri - Columbia, Missouri 65211, USA 

(Dated: November 30, 2009) 

By exploiting our recently derived exact formula for the Lindhard polarization function in the 
presence of Bychkov-Rashba (BR) and Dresselhaus (D) spin-orbit interaction (SOI), we show that 
the interplay of different SOI mechanisms induces highly anisotropic modifications of the static 
dielectric function. We find that under certain circumstances the polarization function exhibits 
doubly-singular behavior, which leads to an intriguing novel phenomenon, beating of Friedel os- 
cillations. This effect is a general feature of systems with BR+D SOI and should be observed in 
structures with a sufficiently strong SOI. 
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Spin-orbit interaction (SOI) is of great interest for 
spintronic applications 0, 0] ■ Electron spin is not con- 
served in the presence of SOI, which allows for purely 
electric manipulation of spins In conjunction with 

other carrier scattering mechanisms, SOI leads to intrigu- 
ing novel phenomena. One of the major findings is the 
detection of the spin Hall effect [f| H, 0, H], predicted 
long ago as an outcome of the interplay between SOI 
and electron- impurity scattering [9]. In turn, electron- 
electron scattering mediates mutual transformations of 
spin and charge currents, occurring due to spin Coulomb 
drag in individual layers [1(3, 11 1 and due to spin Hall 
drag [12] in electronic bilayers. 

In zinc-blende semiconductor nanostructures the inter- 
play between different mechanisms of SOI can itself have 
crucial consequences. In the presence of both Bychkov- 
Rashba (BR) [H] and Dresselhaus (D) Q SOI the sys- 
tem possesses C<z v symmetry. The BR coupling strength 
a depends largely on the asymmetry of structure while 
the D coupling (3 vary mainly with the thickness of struc- 
ture. In the special case when the BR and D SOI 
strengths are adjusted [TBI, [l6| to be equal, even higher 
SU(2) symmetry occurs in the system [17| and various 
relaxation [lj| and optical [l9[ properties of the system 
turned out to be identical to those in the absence of SOI. 
A remarkable demonstration of such suppression of SOI 
is the fresh experimental realization [20| of the persistent 
spin-helix. 

Another distinct manifestation of the interplay of BR 
and D mechanisms is the SOI-induced anisotropy of 
single-particle spectrum, which modifies spin relaxation 
and transport properties of the system 21, H, 23$ . Re- 
cently we have studied the influence of that anisotropy 
on the many-body response of a 2DES 24| ■ Our calcula- 



tions have revealed a fine structure of the plasmon spec- 
trum, which produces a striking asymmetric doublet of 
the structure factor versus momentum orientation. The 
joint action of BR and D SOI leads to dependence of 
the inter- chirality particle-hole continuum on direction. 



Thus, the plasmon propagation may be free in one direc- 
tion, but strongly damped in a different direction, where 
the plasmon dispersion enters the particle-hole contin- 
uum. This creates a possibility of directional plasmon 
filtering, potentially useful for both spintronic and plas- 
monic devices. 

In this Letter we study the many-body response of 
2DES in the presence of BR and D SOI in the static 
limit. By exploiting our recently derived exact formula 
for the Lindhard polarization function (PF), we reveal a 
highly anisotropic behavior of the static dielectric func- 
tion. Particularly, the inter- chirality transitions strongly 
modify the singular behavior of the dielectric function. 
The sharp anisotropic Fermi surface makes the position 
of singularities dependent on momentum orientation and 
on the ratio of the BR and D SOI strengths, in addition 
to the usual dependence on the magnitude of momentum 
q. We find that PF shows a sharp cusp at q < 2kp (kp 
is the Fermi wave vector) for the momentum orientation 
along the [110] direction while in the perpendicular [110] 
direction, the singularity occurs at q > 2k p. Most im- 
portantly, we observe that PF exhibits a doubly-singular 
behavior— the singularities occur both at q < 2kp and 
q > 2k p. As a direct consequence of this, we find that 
the Friedel oscillations propagate with two slightly differ- 
ent spatial frequencies and a novel beating phenomenon 
of Friedel oscillations takes place. At strictly equal SOI 
strengths, a = /3, only the second singularity survives. 
The Friedel oscillations become isotropic and the effect 
of SOI reduces to a simple renormalization of the position 
of singularity. 

The Hamiltonian of BR and D SOI in quantum 
wells of zinc-blende structure, grown on a (001) sur- 
face, is i?soi = a (&xky — o- y k x ) + (3 (a x k x — o y ky) where 
a Xt y are the Pauli matrices, k is the in-plane elec- 
tron momentum with its magnitude k and polar an- 
gle </>k- The eigenvectors of the Hamiltonian H = 
Ha + -Hsoi with H — k 2 /2m* (m* is the electron effec- 
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FIG. 1: (Color Online) (a) The Fermi contour in the presence 
of BR+D SOI. Arrows indicate directions of spin. The thick 
solid and dashed lines are diameters along the [110] direction, 
connecting maximally distant points on the different chirality 
subbands and on the outer subband, respectively. The thick 
dash-dotted line is the inter-subband diameter in [110] direc- 
tion. The BR and D strengths are related as a — 1.21/3 with 
p = 0.2&F. (b) The lengths of respective diameters of the 
Fermi contour as a function of the parameter 9 for p — 0.2kF. 
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beled by the chirality p = ±1. Here A is the nor- 
malization area and the phase of spinor is given by 
ip(a, P, 0k) =Arg[ae^ k + i/3e~^ k ]. The angle-dependent 
BR-D momentum is £(p, 9, 0k) = pyl + sin(2f5) sin(20k) 
where p = to* y 'a 2 + f3 2 . The angle parameter 9 is de- 
fined as tan 9 = [3 /a and describes the relative strength 
of the BR and D SOI. The Fermi momenta of the chiral- 
ity subbands are also angle dependent: k F (p,9, 0k) = 

y/2mEF + £,{p, 9, 0k) 2 — p S,(p, 9, 0k) where the total 
carrier density n determines the Fermi energy, Ep = 
(ttu — p 2 ) /to*. Fig. [1] shows the anisotropic Fermi con- 
tour in the (k Xl k y ) plane. 

In the static limit the dielectric function e(q) = 1 — 
v(q)U(q) where v(q) = 2ne 2 / (noq)F(qd) is the bare 
Coulomb interaction with ko the low frequency dielec- 
tric constant. The form factor F(qd) takes into ac- 
count the transverse width d of the quantum well. It 
goes as 1 — (1/3 — 5/4ir 2 )qd in the long wavelength limit 
qd — ► 0, and as 3/ (4n 2 qd), in the opposite limit qd — > oo. 
The exact PF H(q) can be expressed in terms of the 
non-interacting Lindhard response function, II (q), as 
n((|) = U°(q) [1 - v(q) (1 - G+ (q)) n°(g)r* ■ Here the 
"charge-channel" local field factor G + (q) 25| takes into 
account all electron correlations, related to the vertex 
corrections beyond the random phase approximation. We 
neglect the effect of SOI on G+. 

The static Lindhard PF in the presence of SOI is 
defined as a sum over the indices p and v, H°(q) = 



n° 

fits 



(q) 



dk f F (E^k))-f F (E^(k + q)) 
(2nf 



T v [k,k- 



E„(k) - E„. v {k + q) 

(1) 

Here v — ±1 refers to the intra- and inter-subband 
contributions to H°(q) and fp(E^(k)) are the Fermi 

distribution functions. The form factor T v (jt, k + q\ 

comes from the product of bare vertices and describes 
the overlapping of spinor wave functions. It is given 

by T v (k,k + q\ = \ [1 + ^cos (A<p q )] where we define 

Atfq = ip(a, P, 0k) — <fi(ct, (3, 0k+q)- Notice that in con- 
trast to the case with pure BR or D SOI, here PF depends 
additionally on the momentum orientation, q . Recently 
in Ref. |2J we have derived an exact formula for the PF 
in the presence of BR+D SOI, which is omitted here for 
the sake of brevity. The results, obtained from this for- 
mula are in agreement with the previous classic resul t by 
Stern (2(| and findings by Pletyukhov and Gritsev 27 1, 
respectively, in the limits of vanishing SOI and of pure 
BR SOI. 



Further by exploiting this general formula for PF [24( , 
we calculate the static PF and analyze the modifica- 
tions induced by the BR+D SOI. We use the realistic 
materials parameters for InAs quantum wells by taking 
to* = 0.023mo, ko = 14.55, and the transverse width of 
the quantum well d = 15 nm. We have also defined and 
will use in the following the dimensionless parameters 
x = q/2kp and r = p/kp with kp = \J 2m* Ep + p 2 . 
In Fig. Owe plot II (q) in units of the density of states 
at the Fermi level g — m* /2n as a function of q for 
two different orientations of the momentum, q = 7r/4 
and q = 3"7r/4, and for several values of the angle 
parameter 9. The solid line represents PF in the absence 
of SOI, r = 0, while the dashed line corresponds to PF 
H°brd(\ ( i\) m t ne case 01 t ne pure BR or D SOI, 6 = or 
7r/2, respectively. All the curves of H a (q) show singular- 
behavior at one or two values of q, determined by the 
lengths of diameters in Fig. QJi and [TJd. Independent of 
9, the curves, which refer to q = 37r/4, exhibit a singu- 
larity at the wave vectors q\ c = 2kp\ / l — r 2 sin 29 < 2kp 
(cf. the solid line diameter in Fig. [T|) with a max- 
imum polarizability at the singular point, LT 1 nax , 
exceeding the maximum value of n'g^dcfl). At 
wave vectors qi c > 2kp PF develops a second point 
of nonanalyticity with a local maximum n™ ax < 2. 
Here q2 C = niax0 k Q(0k) with the function Q(0k) = 

(sin(0 k ) - cos(0 k )) £ k + V 1 + r2 sin ( 2( M sin (29) 
(cf. the dashed line diameter in Fig. [lj. For values 
of 9 not far from 7r/4, the singularity at qi c becomes 
well pronounced while the singularity at q\ c turns into 
a sharp cusp. The curves in Fig. [21 which refer to 
<p q = 7r/4, exhibit only one singularity at the position 
<73 C = 2kp\J\ + r 2 sin 29 > 2kp (cf. the dash-dotted line 



3 



2.0 






no SOI iji" 

. 6 — 0,n/2 \V, 




V {a) 




(9 — 7r/16 \ 




1.8 


S = n/S ^ 

6 = 7?r/32 


v. 




■ «= 15^/64 


\Vv 


1.6 


fl = ir/4 

i . i 










0.98 1.00 


1.02 




ql2k F 






FIG. 2: (Color online) The static polarization function II (q) 
as a function of the momentum q for its two orientations: (a) 
0q = 37r/4 and (b) 7r/4. The different curves correspond to 
different values of 9, shown on the graph legends for r = 0.1. 



diameter in Fig. [T]) and the polarizability in the singular 
point varies within the 2 < U°(q) < T1'bhd(\'i\) window. 
In all singular points PF is continuous and its derivative 
discontinuous. At the special values of 9 = tt/4 and 
37r/4 when a = ±/3, the form factor T v ceases to depend 
on the wave vector k and reduces to the Kronecker 
symbol. In this case we find that the effect of SOI 
on the density response of a 2DES disappears for any 
value of (j)q. The only remaining modification reduces 
to a simple renormalization of the isotropic Fermi wave 
vector, kp — > fcj 



Fc = VT+ r 2 k F . 
The electron density deviation, generated by the per- 
turbation of a single impurity, which is embedded at 
R = in the electron sheet in the (x, y) plane, is de- 
termined by the static density response function via the 
relation 



Sn{R) 



dq 
(2^ 



(2) 



where Vi(q) is the Fourier transform of an isotropic im- 
purity potential. Since in the presence of BR+D SOI PF 
is anisotropic, we can rewrite Eq. [2] as 



5n(R) 



dxxVi (x) 
o Jo 



^i£x cos( 



^n(x,0 q ) 

where £ = 2kpR = (2^/2/r s ) (R/as) and </>r is the 
polar angle of the vector R, r s — \/2/kpaB the di- 
mensionless constant of electron-electron interaction [25| , 
as the effective Bohr radius. At large distances from 
the impurity, the integrand has a rapidly-varying phase, 
£x cos (</) q — </>r) so the main contribution to the inte- 
gral comes from the point where the phase is station- 
ary and the exponential function e l t x cos (<P<i- <Pr) oscillates 
less rapidly. The application of the method of stationary 
phase to the integral over </> q in ([3]) yields 



Sn(R) 



4: \f^ J dxy/xV imp (x)TL(x, 0r) cos - £x 

(4) 



Here we have used the relation H(x, 0) = H(x, 7r), which 
is the case in the presence of BR+D SOI due to C2 V 
symmetry. Thus, at large distances from the impurity 
one can replace the second argument of PF </> q by (/>r. 

In order to take the integration over x we can exploit 
the Riemann-Lebesgue lemma, which says that if a func- 
tion oscillates rapidly around zero then the integral of 
this function is small and the principal contribution to 
the integral is determined by the integrand behavior in 
the neighborhood of singular points. Therefore, one can 
reduce Eq. [5] to the following asymptotic expression for 
the density deviation 



8n{R) 



where 



(5) 



^9 Vimp i&ic ) 



[1 - v(x lc ) (i - G+ (x ic )) n°{x lc , <f> R )y 



(6) 



and Xi C = Xi C (r, 0, q ) denotes the position of the ith 
singularity of PF. The increment of PF SU°(x, cc,- c ) = 
H°(x,4>ji) — n°(xi C ,0R) near the singularity be 
represented as 



5H°(x, xl 



-g-d [± (x 



2 \ a ±,i 



(7) 



where i9 (x) is the unit step function, the signs ± cor- 
responds to the function below (x < Xi C ) and above 
(x > Xic) the singularity Xi c . The critical exponents a± t i 
and the coefficients a i describe the power law behavior 
and the maximum polarizability at the singular points. 
In Eq. [5] we have assumed that the nonanalytic behav- 
ior of the interacting PF Tl(q) is completely determined 
by its noninteracting part H°{q) [3]. Substituting Eq. [7J 
into Eq. [5] and making use the Lighthill theorem [281 ] . 
after the integration over x we get 



cos 



Sn(R) 2A ^ (a± ,)! 

V * ±,l 

(8) 

As we have already discussed, PF can qualitatively 
change its singular behavior, depending on the ratio of 
the BR and D SOI coupling strengths as well as on the 
momentum orientation. Accordingly, the Friedel oscilla- 
tions, given by Eq. [8l can exhibit completely new fea- 
tures. For instance, at = ir/8 PF for two orthogonal 
orientations, q = 7r/4 and </> q = 37r/4, shows singu- 
larities, respectively, at the positions x+ c and x- c with 
x± c ~ 1 ± 5c- It is clear that the phase difference be- 
tween these two orientations is about (x +c ~ X—c) £ — 
(4v / 2/r s ) S c {R/as) and can result in a striking differ- 
ence in the behavior of Friedel oscillations at the distance 
R/as of the order of r s /A^/28 c - In InAs samples with the 
electron density n = 10 16 m~ 2 , we have r s ss 0.12 and for 
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FIG. 3: (a) The doubly-singular polarization function for a = 
1.21/3 with p = 0.2&F (solid line). Inset shows the height of 
cusp vs 8. The dotted and dashed curves represent II (q) for 
a = f3 with p — 0.2fci? and p = 0, respectively, (b) Beating of 
Friedel oscillations, induced by the doubly-singular behavior 
of n°(g) in (a). The electron density n = 2 x 10 16 m" 2 . 



r = 0.1 taking S c = 0.5 x 10~ 2 we obtain that the Friedel 
oscillations in the <^>r, = 7r/4 and </>r = 37r/4 directions 
are in antiphase at the distance of the order of R < 5a b- 
Another interesting effect appears when PF exhibits 
the doubly-singular behavior (cf. the curves with 9 = 
7?r/32 or 7?r/32 in Fig. H}. In Fig. [3h. we illustrate the 
doubly-singular behavior of PF separately for r = 0.2. 
Inset shows the height of the cusp at q = qi c (cf. the 
solid curve in Fig. Q})) versus 9. As seen, n5 nax = 11^^^, 
at 9 = and increases with 9 up to its maximum 
ITf ax « 2.17 at about 9 = 0.18tt. With a further increase 
of 9, n™ ax drops sharply to its value in the absence of 
SOI, ITf ax = 2 at 9 = vr/4. By fitting numerically the sin- 
gularities of the exact PF, we have established with nu- 
merical certainty that the critical exponents in the power 
law behavior in Eq. [7] are (cf. Ref. 1271): ot-,i,i = 1 and 
ct+,1,2 = 1/2, respectively, for the left- and right-side of 
both singularities at x\ c = 1 — 5i c and X2 C = 1 + #2c- 
Hence, in the limit of large £ the contributions to Sn(R), 
coming from the left-side neighborhood of both singulari- 
ties, are small. Taking into account also that the smooth 
functions A(xi C ) and a i do not depend on the large pa- 
rameter £, one can approximate a i A(xi C ) ~ 2A(1) = Aq. 
Thus, the density deviation at large distances from the 
impurity can be reduced to the following simple form 



5n(R) 



-£2 l sm 0- 



(5i c K + sin(l + <5i c ) f] 



(9) 



It is clear that due to the existence of two singularities 
at x\c = 1 — Sic and X2 C = 1 + #2c, the Friedel oscilla- 
tions sin (1 — Sic) £ and sin (1 + 5 lc ) £ propagate with two 
quite close spatial frequencies and a beating phenomenon 
of Friedel oscillations can be observed at the beat fre- 
quency \/2 (Sic + S2c) /r s . Fig. [3Jd illustrates the first 
two destructive interferences of the Friedel oscillations 
that occur at R ~ 3a b and 9a b- Taking into account 
higher order terms in £ and p will partially smooth the 
interference picture, however, the beating of Friedel oscil- 
lations as a distinct modulation of the density deviation 



is a stable feature of systems with BR+D SOI and should 
be observable in experiment. Notice also that in samples 
with a stronger SOI such as HgTe quantum wells, the 
separation 5\ c + S2 C between the singularities increases, 
which will essentially facilitate the experimental detec- 
tion of the destructive interferences. 

In conclusion, we have calculated the static response 
of a 2DES in the presence of joint BR and D SOI. We 
find that one of the main modifications is the induced 
shift of the singularity position of the static PF, which 
is in opposite directions for orthogonal momentum 
orientations. This results in a strong anisotropy of the 
Friedel oscillations. More interestingly, we have shown 
that in certain situations PF exhibits a doubly-singular 
behavior, which generates a novel phenomenon - the 
beating of Friedel oscillations. This intriguing prediction 
exemplifies how usually weak SOI can generate a 
qualitatively new and physically robust occurrence as a 
measurable signature of the many-body response of a 
2DES. 
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